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ABSTRACT 

Bimetric MOND (BIMOND) theories, propounded recently, predict peculiar gravi- 
tational interactions between matter and twin matter (TM). Twin matter is the hypo- 
thetical matter that might couple directly only to the second metric of the theory, as 
standard matter couples only to the first. Considerations of cosmology in the BIMOND 
framework suggest that such TM might exist and copy matter in its attributes. Here I 
investigate the indirect interactions that BIMOND theories predict between local, non- 
relativistic mass concentrations of matter and TM. The most salient result is that in 
the deep-MOND regime of the matter-TM-symmetric theories, TM behaves as if it has 
a negative gravitational mass relative to matter (active and passive, with the inertial 
mass still positive). To wit, interaction within each sector is attractive MOND gravity, 
but between matter and TM it is repulsive MOND gravity. Using the space-conformal 
invariance of the theory in the deep-MOND regime, I derive various exact results; e.g., 
the repulsive force between a matter and TM point masses (space-conformal theories 
are a natural framework for masses of opposite signs). In the high-acceleration regime, 
the interaction depends on a parameter, j3 (the strength of the Einstein-Hilbert action 
for matter). For the favored value f3 = 1, matter and TM do not interact in this regime; 
for j3 < 1 they attract; and for f3 > 1 they repel each other. Such interactions may have 
substantial ramifications for all aspects of structure formation, such as matter distribu- 
tion, peculiar velocities, and effects on the CMB. The repulsive interactions probably 
lead to segregation of matter and TM structures, leading, in turn, to intermeshing of 
the respective cosmic webs, with high-density nodes of one sector residing in the voids of 
the other (possibly conducing to efficient evacuation of the voids). Weak gravitational 
lensing by TM seems the best way to detect it or constrain its attributes. In the MOND 
regime a TM body acts on matter photons as a diverging lens. Strong lensing occurs in 
the high acceleration regime, and thus depends on f3. For f3 = 1, a TM mass does not 
bend (matter) light in the high-acceleration regime: no strong lensing effects of TM are 
expected in this case. I also discuss briefly asymmetric theories. 



Subject headings: 
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1. Introduction 

Bimetric theories of gravity involve two metrics as independent degrees of freedom: one felt 
directly by standard matter, g^, and an auxiliary metric, g^ v . Aspects of such theories have been 
extensively discusses; e.g., by Isham, Salam, & Strathdee (1971), by Rosen (1974), and, for some 
recent treatments, with references to other and to earlier work, see, e.g., Boulanger & al. (2001), 
Damour &; Kogan (2002), Bias, Deffayet, & Garriga (2006), Bahados, Ferreira, &: Skordis (2009), 
and Bahados, & al. (2009). It may be the case that the metric g^ is indeed merely an auxiliary 
field in the description of gravity of standard matter. It is, however, natural to suppose that g^ v 
comes with a matter sector of its own, and should be viewed on a par with g^ v . This opens the 
way for introducing twin matter (TM), which may couple to its ilk in the same way as standard 
matter does (electromagnetically, weakly, etc.), but which couples gravitationally only to <? M „, just 
as matter couples only to g^ v . Since the two metrics are coupled through a gravitational term in 
the action, matter and TM do couple indirectly by some unconventional gravitational interaction. 

A bimetric formulation of MOND (BIMOND) has been propounded recently (Milgrom 2009b), 
which points even more forcibly to the possible existence of TM. 

BIMOND is governed by an action of the form 
I = ~Y^q I W 2 R + ag^R- 2{gg) 1 ' A f( K )a 2 M(r / a^x + I^g^,^) + I M (g, v , X i)- (1) 

Here T represents a collection of scalar variables formed by contractions of the acceleration-like 
tensors C^ 7 = Tp—Tp, where and r^ 7 are the Levi-Civita connections of the two metrics. Also, 
k = (g/g) 1 ^ (g and g are minus the determinants of the two metrics), G is the phenomenological 
Newton constant, and I use unites where c = 1. (Ai may also depend on scalars constructed 
from the two metrics, such as k.) I m is the action for matter, whose degrees of freedom and their 
derivatives are collectively marked by ipf, these interact among themselves, and couple to the metric 
g^ v alone in the standard way, and similarly I M is the TM action. I have worked, in particular, 
with favorite choices of a the scalar variable in M built from the tensor 

for example, 

T = —gTT^. (3) 

To obviate possible confusion, especially in the present MOND context, I note at the outset 
that the TM is not the putative dark matter (DM), and is not taken to play its role of enhancing 
gravity in galactic systems. It is still MOND departure from standard general relativity (GR) that 
replaces dark matter. TM may still produce some effects, e.g., in structure formation, that are 
conventionally attributed to cosmological DM. It may also linger in otherwise matter-dominated 
territory to produce some visible effects (see 4.0.3). 
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We have no a priori idea of how the TM sector looks like; neither do we have any observational 
information on the subject. Does it exist at all? Is it made of the same 'stuff' as matter, and is it 
subject to the same physics? Is it present in the Universe in the same amounts? Has it undergone 
similar processes in cosmic history (big bang, inflation, seed fluctuations, structure formation, etc.)? 
We do not even know whether matter coupling to gravity is the same in the two sectors. To be able 
to progress I shall assume symmetry between the sectors to the effect that TM duplicates matter 
in its composition, interactions, etc., and that if (g^,Xi) is the same configuration as (g^,tpi), 
then 1^1(9^,^1) = {oi/ j3)I M {% v ^ Xi)- This ensures that in the absence of the interaction Ai, the 
dynamics within the two sectors are identical. With the interaction the above assumptions do not 
ensure complete symmetry when «//3. 

In cosmology, the above assumptions plus the assumptions that the cosmic matter contents 
of the two sectors are the same, and a symmetric choice of f(n) = ((3k + an" 1 ) /(a + /3), was 
shown (Milgrom 2009b) to lead to cosmologies in which g^ v = g^, with each metric describing 
a standard, Friedmann- Robertson- Walker (FRW) cosmology, with a cosmological constant A = 
—a^M.(Q)/(a+j3). This automatically retains the known successes of such a cosmology. Cosmology 
is thus, arguably, the strongest motivation for postulating that TM exists. In such BIMOND 
cosmologies, G//3 appears in the matter sector of such a cosmology as Newton's constant. This 
may induce us to prefer ^wl. The value f3 = 1 is also special because such BIMOND theories have 
a simple limit when a — > 0: they tend to GR in the matter sector when M'(z) — > for z — > 00, 
which is required for the nonrelativistic (NR) limit to yield Newtonian dynamics for a — > 0. I have 
not yet investigated the limit a — > for BIMOND theories with /3 / 1 to see to what extent they 
differ from GR. I will, none the less, keep the discussion more general, taking /3 = 1, as an example, 
only at a later point. 

While on cosmological average the two metrics might be the same, departures from equality 
must occur due to random fluctuations in density, which are different in the two sectors. To treat 
large-scale-structure formation through the development of small perturbations on the background 
of the expanding Universe (including their imprints on the CMB, etc.) we need to expand the 
BIMOND equations of motion around the double FRW metric describing cosmology at large. 

Note, importantly in this context, that even if early on the Universe is characterized by high 
accelerations (as measured, e.g., by cH, where H is the Hubble parameter), what determines 
whether we are in the MOND regime is the argument of the interaction function M.. This argument 
is small when the two metrics are near each other. Thus, structure formation in the initial stages, 
in an otherwise symmetric universe, occurs in the MOND regime. So, even at small amplitude of 
fluctuations the development is nonlinear, as the MOND potentials are not linear in the density 
fluctuations, right from the outset. Structure formation, in all its aspects, such as the distribution 
of matter, peculiar motions, and CMB fluctuations, are thus phenomena that might be greatly 
affected by the new physics inherent in BIMOND, in particular by the possible effects of TM. The 
treatment of this problem has to be left to numerical simulations. 
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In this paper I take up the more modest task of investigating the dynamics of well formed 
systems involving well separated masses of matter and TM, of sizes much smaller than cosmic scales. 
These are then taken to be NR systems on a double Minkowski background (<? M „ ~ g^ v ~ rj^)- The 
appropriate limit of BIMOND theories under these circumstances was shown to be the class of NR 
MOND theories described in detail in Milgrom (2009c). 

The analysis here includes derivation of the forces between combinations of matter and TM 
bodies, and also the characteristics of gravitational lensing of matter photons by TM. 

In section 2, I describe the field equations that govern NR interactions within and between 
the two matter sectors. In section 3, I specialize to the class of theories that are fully symmetric in 
matter and TM. In section 4, I consider some possible observational consequences of the existence 
of TM and its peculiar interaction with matter, including gravitational lensing. In section 5, I look 
briefly at examples of asymmetric theories. Section 6 is a discussion. 



2. Formalism 

For a system of slowly moving matter and TM distributions p and p, respectively, on a double 
Minkowski background, for the choice of variable Y as in eq.(3), the solution for the metrics of the 
BIMOND field equations is (in some gauge) 

fiW = ~ 20<5 M „, <? M „ = f]^ - 2# M „ (4) 

(Milgrom 2009b). 1 Here, the potentials <fr and <p are solutions of the NR limit of the BIMOND 
equations; these equations can be derived from the Lagrangian L = J Cd 3 r, with 

C = ~&b {a( ^ )2 + ~ °% M ^<I> ~ W)V«o]} + P(^v 2 - 0) + 1P {U 2 - 4>). (5) 

The last two terms in the Lagrangian density tell us that <fi is the (MOND) gravitational potential 
for matter, and cf> is that for TM, in the sense that the acceleration of matter and TM test particles 
is given by v = — V</>, and v = — V0, respectively. 2 

Invariance of the theory to translations implies that for a closed system we have a vanishing 
total force: 

F = - / {pVcp + 7/ 3V^)d 3 r = 0, (6) 



1 Choices of the acceleration scalar arguments of M other than that given in eq.(3) result in other NR limits of 
the corresponding BIMOND theory than those given by eq.(4). In general, more potentials are needed to describe 
each metric, and a more complicated set of coupled equations for these have to solved (Milgrom 2009b). 

2 This NR theory is, in fact, the NR limit of a group of variations on the relativistic theory, where we can use in 
the contraction in T, g**" instead of g 1 *" , or use combinations of such variables, for more symmetry between the two 
metrics. 
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and the conserved momentum is 



J (pv + jpv)d 3 r. (7) 



I define p as nonnegative, so the sign of 7 matters. In the spirit of what I said above about 
the relativistic theory, I take from now on 7 = a//3. 3 

The field equations are 

A<j) = 47rG/rV + /r x V • (M'V0*) = 4ttG(p + Pp ), 

A<f> = AirGP^p - a^V ■ (M'V(f>*) = 4nG(p + p p ), (8) 

where 4>* = 4> — (f>, and p p and p p playing the role of 'phantom matter' (PM) densities for the two 
sectors, (a = 0, or j3 = do not give MOND theories, and I exclude such values.) Note that for a 
given configuration the amounts of PM felt by matter and TM are, in general, different. 

Subtracting one equation from the other we get a decoupled equation for 4>* 

V>*(|V^/a )V<f]=47rG(p-p), (9) 

where 

p*(x) = P-^±1m\x 2 ). (10) 
a 

There is an equivalent, but more transparent, way to write the theory for the case a + j3 7^ 0: 
Define 

M(z) = -M(z/q) + T ^-, (11) 
where q = a /(a + (3) 2 . Then, the La grangian density can be written in terms of the potentials 

^ = /^ + a 0, = aC(0-0), (12) 

where ( = (a + /3) _1 , as 

C = ~8^G {C( ^ )2 + a oM(^)V«o]} " C<Xp + - Kp -p) + \{pv 2 + y^ 2 ). (13) 

The two potentials are now decoupled, satisfying the equations 

A0 = 4vrG(p + ^p), V • {M'[(V4>/a o ) 2 ]V0} = 4vrG(p - p). (14) 



The potential 4> is thus a linear combination 



= r + (15) 



3 In the present context, we can absorb \a//3\ in p; so, this choice of I7I may be viewed as the choice of a convenient 
normalization for p. 
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where the Newtonian potentials in the two sectors, cf> N and <f> N , are the solutions of the Poisson 
equation for p and p separately. The matter and TM MOND potentials are then gotten as linear 
combinations 

<f> = (4>+4> = (l+\)- 1 (a- 1 cp N +f3- 1 N )+0, j> = C$--4> = (l+X)- 1 (a' 1 cp N +(3- 1 ^ N )-^4>, (16) 

a 

where A = f3/a. 

The deep-MOND limit is formally implemented by taking a — > oo, G — > 0, with Ga kept 
fixed. In this limit, 4> becomes super-dominant over the Newtonian-like potential (p. We then have 
to have M'(w) — > w 1/2 giving for the matter potential 

V-(|V0|V<A)=47rGa o (p- / 5), (17) 

and 4> = —(f3/a)4> for the TM potential, so we can also write 

V • (|W|V0) = 4irGa (p-p). (18) 

If we put matter and TM test particles at the same position, r, in a gravitational field, the 
former will have an acceleration a = — V0(r), while the latter will be accelerated by —((3/a)a. 
(Matter and TM do not follow the same weak equivalence principle, as is clear from the very 
construction of BIMOND, with matter following geodesies of g M „ and TM those of g^„ ) In the 
deep-MOND limit of theories with a and (3 of the same sign, matter attracts matter as in standard 
MOND; TM attracts TM as in MOND, but with an effective value of Ga Ga (/3/a) 2 ; however, 
matter and TM repel each other. 

To get the Newtonian behavior in the matter sector for a — > 0; namely, to get (ft — > (j) N in this 
limit, for a pure-matter system, 4 we have to have 5 •M'(oo) = 1/(1 — (). 

Since Ai'(w) « w 1/2 for small w > 0, and is thus necessarily positive there, we cannot have 
C > 1, lest M'(w) vanishes at a finite value w, which it must not. This is the condition derived in 
Milgrom (2009c) on the same grounds. We can have C = 1 (and then M' — > oo at infinity); but 
this is excluded by solar-system constraints. There may be other constraints on the values of a, f3 
from different consistency requirements in the relativistic and NR theories, but such are yet to be 
found. 

With the above value of Ai'(oo), we have in the high-acceleration regime 



4 Note that the Newtonian limit is defined by the requirement that it reproduces Newtonian dynamics in the 
Matter sector, not necessarily in the TM sector. 

5 This has to be required only because we insisted that G is Newton's constant. Had we started from some general 
coupling, G' , this relation would just constitute a definition of G in terms of G' . 



1 + 



a/3 



a 



V 



(19) 
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We see that while G is the gravitational constant in the Newtonian limit of the matter sector, 
it is G = G[l + (ft — l)(ft — a)/a(5] that plays this role in the TM sector. For ft = 1 or ft = a we 
have G = G but otherwise they are different. In fact, for some choices of a, ft we have G < 0, in 
which case the a — >■ limit corresponds to repulsive gravity in the TM sector. If we deem this 
undesirable we can eliminate the corresponding a, ft values. (For a, ft > we have G > with our 
already assumed inequality ( < 1.) 

We also see from eq.(19) that for ft = 1, (f> = 4> N and (f> = cf) N , which means that there is no 
interaction between the two sectors in the Newtonian regime. This is true for the fully relativistic 
theory, where for ft = 1 (irrespective of a) BIMOND separates in the limit a ( , — > to two copies of 
GR in the two sectors [possibly with a cosmological constant ~ a%M.(oo)]. 

When a ^ ft the dynamics in the two sectors can be quite different. Such theories are worth 
investigating, but my purpose here is not to conduct an exhaustive study of this class of theories, 
only to demonstrate some salient results. To this end I shall consider mainly theories that are fully 
symmetric in matter-TM, namely those with a = ft. I will then discuss briefly some asymmetric 
cases. 



2.1. A quasi-linear formulation 

Solution of the field equations (14) requires solving a nonlinear Poisson equation for the given 
matter-TM configuration at hand. This may be rather taxing, especially when applying the theory 
to time-dependent problems, such as that of large-scale-structure formation. As in the case of 
the quasi-linear formulation of MOND (QUMOND; Milgrom 2009c), which parallels the Modified- 
Poisson-equation formulation of Bekenstein & Milgrom (1984), I describe here a quasi-linear theory, 
derivable from an action, not equivalent to our theory, but which captures much of its essence, and 
which should be much easier to apply. As in the case of QUMOND, this requires adding an auxiliary 
potential ijj to the MOND potentials for matter and TM. 6 Consider the Lagrangian density 

C = -^{ 2 «C(W-v#vv>+C(/3V^ (20) 

In terms of the two potentials (f> and (p related to the MOND potentials by eq.(12) one can write 

C = ~8^G {2 ^^ ' ^ + C( ^ )2 " a « Q t( W°») 2 ]} " -P)~ + + \p w2 + Yft^' (21) 
giving the field equations 

A4> = 4irG(p + ^p), AiP = 4irG(p- p), = V • {Q'[(VV/a ) 2 ]VV}. (22) 



6 Since (f> and <f> are decoupled in the Lagrangian density in eq.(13), we simply apply to <f> the same procedure that 
has lead to QUMOND. 
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It is easy to see that with the appropriate choice of Q [namely, with Ai'(x 2 )x = y being equivalent 
to Q'(y 2 )y = x] the theory is a very good mimic of equations (14). For example, in a spherically 
symmetric case they are identical. 

3. Dynamics in the fully symmetric theory 

In the fully matter-TM symmetric case, a = ft, which seems to give appealing cosmological 
solutions, we can write the field equations (14) as 

A0 = 4irG(p + p), V ■{M'[{V4>/a ) 2 )V4>} = ^G{p- p), (23) 

and the MOND potentials are then given by 

= C0 + <A, = C0-<A, (24) 

where now ( = (2ft)" 1 . 

In the deep-MOND limit <fi satisfies eq.(17), and 4> = — <fi. 

Interestingly, in the deep-MOND limit of the theory, TM behaves as if it has a negative active 
and passive gravitational mass relative to matter, while its inertial mass is still positive: Negative 
active mass because p enters the source for the matter 4> potential with a negative sign, and negative 
passive mass, because it is accelerated by the gravitational field <p = —(f). So bodies in the same 
sector attract each other, while bodies in different sectors repel each other. I discuss gravity in this 
important MOND limit in more detail in subsection 3.2. 

The case p ~ p has to be commented on. In this case (f> ~ 0, while <p ~ 2<j) N ~ 2<fi N . The 
MOND limit applies when [Ga (p — p)R] 1 ^ 2 S> G(p + p)R, where R is the characteristic size of 
the system. The configuration of density near-equality is unstable, and with a small separation, 
repulsion occurs (see subsection 3. 3). 7 

Consider now the Newtonian limit, a — > 0, where from the general eq.(19) we have 

= + (2(-l)4> N , 4> = 4> N + (2(-l)r. (25) 

We see that, unlike the MOND limit, here the predicted fields do depend on (. Recall that 
for the choice, £ = 1/2 (ft = 1), we have (f) = (f> N , <p = <f> N , with each sector seeing exactly its own 



7 A similar result applies in the more general, relativistic case: In configurations where the energy-momentum 
tensors in the two sectors are equal (with our normalization 7 = a//3) the solution of the field equations is p M „ = g^, 
with both being the solution of the standard Einstein equation for the configuration [with a cosmological constant 
oc doM(0)]. For example, a double Schwarzschild metric is a spherically symmetric, vacuum solution of the BIMOND 
equations. Matching it to an interior solution will show that it corresponds to a central mass made of equal amounts 
of matter and TM. 
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Newtonian potential; so, matter and TM do not interact at all. For £ > 1/2 matter sees its own 
Newtonian potential plus a fraction 2£ — 1 > of that of TM; so they attract each other with a 
reduced effective gravitational constant. For ( < 1/2 they repel each other. 

The case £ — > (/? — > oo) is also interesting: We see from eq.(24) that in this case (p = ~ ( t ) = 
<f) for the full range of the theory (Newtonian-MOND); <f> becomes immaterial, and the MOND 
potential <p satisfies the second of eq.(23). We saw that for Q = we have A-f (oo) = 1. So in this 
case we end up with the MOND theory of Bekenstein & Milgrom (1984) with p,(x) = M'(x 2 ), but 
with TM entering the theory as having a negative gravitational mass. 

3.1. Forces on bodies 

The force F (F) on a matter (TM) body that constitutes a subsystem of the density p (p) 
within the volume v (v) is 

F = — / pV(f>d 3 r, F = - / P V4>d 3 r. (26) 

Because the theory is nonlinear we cannot use in these expressions the potential produced by the 
system with the body in question excluded (as can be done in the linear case). As a result, even 
for a point mass m at position r we cannot write the force simply as F = — mV^(r), where (f> is 
produced by the rest of the system. In fact, the force is not even linear in the mass of the body, 
and becomes so only for test particles. In Milgrom (1997, 2002a), I discussed in detail general 
properties of forces in such nonlinear theories. Here I essentially use the results from these papers. 

On dimensional grounds we can write the force between two matter point masses, M and m, 
a distance r apart, as F = —a Mf(m/M,r/R M ), where R M = (MG /a ) 1 / 2 is the MOND radius 
for the mass M, and a negative sign signifies attraction. The same expression applies to the force 
between two TM point masses. For a matter point mass M and a TM point mass M we can write 
F = a Mf*(M/M,r/R M ). 

In the test-particle limit, where one of the masses is much smaller than the other, f(q, A) 
and f*(q,X) are easy to obtain for all values of A, because the force on a test particle equals its 
mass times the gradient of the potential produced by the massive ponit mass, which can be gotten 
analytically. In the Newtonian limit A < 1 we see form eq.(25) that f(q, A <C 1) ~ qX 2 , which 
reproduces the Newtonian force, and f*(q, A ^ 1) « (2( — l)q\~ 2 . Expressions for these functions 
in the deep-MOND limit A> 1 will be given in subsection 3.2 

3.2. Gravitational interactions in the deep-MOND regime 

I now discuss in more detail matter-TM gravity in the deep-MOND regime of the symmetric 
theories. Look at a system of matter (p) and TM (p) where the surface densities are everywhere 
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small enough that the accelerations everywhere are much smaller than a . In this case, the potential 
field is determined from eq.(17), and the potential felt by matter is 4> itself while that felt by TM is 
—(f). To the MOND forces determined from these potentials one adds the subordinate Newtonian 
forces, coming from the (j) potential, weighted by £, as per eq.(24). 

Consider now forces on non-test bodies of both types. In Milgrom (1997), I showed that 
the theory described by eq.(17) is invariant under conformal transformations, just as the two- 
dimensional Poisson equation is. This enables us to derive some very useful exact results, which are 
otherwise difficult to obtain in nonlinear theories of the type discussed here. I also note, in passing, 
that the full symmetry of the theory is then the same as the isometry of a de Sitter space-time, 
with possible ramifications discussed in Milgrom (2009d). In addition to the usual invariance to 
translations and rotations, which have the usual consequences (e.g. conservation of momentum and 
angular momentum) , the conformal symmetry in 3-dimensional Euclidean space includes inversions 
about a sphere of any radius R, centered at any point r . Namely, transformations of the form 
r — > r' = r„ + R 2 \r — r„|~ 2 (r — r ); and also dilatations r — > r' = Ar. The symmetry means that 
if <fi(r) is the MOND potential for mass distributions p(r), p(r), then </>*(r') = <fi(r) is the MOND 
potential for mass distributions p'(r') = J~ 1 p(r), p*(r') = J _1 /5(r), where J(r') is the Jacobian 
of the transformation J = \\dr' /dr\\. For example, from dilatation invariance follows that if 0(r) 
is the MOND potential for p(r), p(r), then 0(r/A) is the MOND potential for mass distributions 
A- 3 p(r/A), A- 3 p(r/A). 

The existence of gravitational masses with opposite signs fits very naturally in such a confor- 
mally invariant theory. To take advantage of the symmetry, the point at infinity has to be treated 
on a par with all other points, since inversions, which are conformal transformations, interchange 
the point at infinity with a finite point. This implies that we should view the Euclidean space as 
the 3-dimensional sphere from the topological point of view. On the other hand, the gravitational 
field lines going to infinity away from a system of finite total mass, M, converge at infinity with a 
net flux, implying an effective negative mass —M, there. Under inversions this mass is brought to a 
finite point, leaving us with a configuration involving a negative mass, even if we start without one. 
So, having such masses from the start is an advantage. Put differently: a compact manifold such 
as the Euclidean sphere, as our space has to be viewed in a conformally invariant theory (here, to 
space transformation, not space-time), cannot accommodate a finite total mass: Applying Gauss 
theorem to any surface that separates space in two, implies that the total mass on one side has to 
be equal and opposite that on the other side, hence negative masses are needed. 

Next, I describe several corollaries of the conformal symmetry. They pertain to a system made 
of matter and TM point masses, m„ m k (all defined as positive), at positions r u r fc , the forces on 



8 For our results to apply it is enough that the accelerations are small across the typical inter-particle distance in 
the system. For example the bodies we treat may be point masses, so the acceleration field near them is higher than 
a a . All we require then is that the Newtonian regimes of the different masses stay far away from each other. The 
point masses can then be taken as spheres larger then their MOND radii. 



- 11 - 



which are F i5 F k . Translation invariance dictates ^\ F 4 + ^ fc F fc = 0, and rotational invariance 
implies r.xF.I £) fc r fc x F k = 0. 

(i) A 'virial relation' holds, which takes the form 

J> • F 4 + J> • F fc = -\(a a G)^[\ ^m, - J> fc | 3 / 2 - 5> 3 / 2 - ^m 3 / 2 ]. (27) 

There is an additional relation that applies for systems with vanishing total 'charge' (see Milgrom 
1997). 

(ii) The following are corollaries of relation (27) for the two body case: For two matter masses, 
say M at the origin and m at r, the force on m is 

F = -|(a G) 1/2 [(M + mf' 2 - M 3/2 - m 3/2 }^; (28) 

so /(g, A>l)w (2/3)A _1 [(1 + q) 3/2 — 1 — q 3/2 ]. This force, which is the same for two TM masse, is 
always attractive. For a matter mass M at the origin and a TM mass M at r, the (repulsive) force 
on M is 

F = -(a G) 1/2 [M 3/2 + M 3 ' 2 - |M - M| 3/2 ]4; (29) 

o V 

so, f*(q, A > 1) » (2/3)A" 1 [l + 9 3/2 - (1 -g) 3/2 ]. For M = M, F = (4/3)(a„G) 1/2 M 3 / 2 (r/r 2 ). This 
is stronger by a factor 1/(2 1/2 — 1) « 2.4 than the attracting force between two equal matter or TM 
masses. 

(iii) The theory is nonlinear; so forces are not additive: When we have more bodies present, 
the force on each has to be determined through a new calculation of the whole system. This, to my 
knowledge, is impossible to derive analytically, in general, but there are exception, several examples 
are discussed below to demonstrate that matter-TM repulsion is general (see also Milgrom 2002a 
for general theorems on this issue). 

(iv) The forces in a three-body system with 'vanishing total charge'; i.e., with ^ m ; = ^ fe m k , 
was calculated in Milgrom (1997). The force on one of the bodies, call it 1, is 

fx = /3i2 t— ^— — — 72 + ftaA-^ (3°) 
\r 2 -r 1 \ 2 |r 3 — r x | 2 

where (3^ = (2/3)(a G , ) 1//2 (|g , j+g , i | 3/2 — |(j , i| 3/2 — |<Zj| 3//2 ) ) with q i = m ; for a matter body, and q k = —fn k 
for a TM body. So, interestingly, in this case the force is simply the sum of the two forces that 
would have applied by the two bodies separately. 

(v) Other configurations for which there are analytic results involve symmetric configuration. 
For example, take a configuration of ./V equal masses M, of the same type, symmetrically placed at 
a distance r from the origin, so that the forces on all are equal and radial (such as at the corners 
of a square, a cube, etc., or forming a uniform spherical shell in the limit of large N). Application 
of relation (27) give the force on each mass as — (2/3)M(MGa ) 1 / 2 (N 3 / 2 — N)rr~ 2 . Now place any 
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mass m of the opposite type at the center and the force becomes —(2/3)(Ga ) 1 / 2 (\NM — m| 3//2 — 
— m 3/2 )rr 2 . The addition of m clearly weakens the self attraction towards the origin, and 
reverses the sign of the force for large enough m. 

In a similar fashion, if we have an equal number of matter and TM point bodies of the same 
mass arranged symmetrically, so that all forces act radially and are equal (for example two masses 
of each type on alternate corners of a square), relation (27) tells us that they are all repelled from 
the center. 

(vi) The potential field of a pair of equal matter and TM masses, M, at r 1 and r 2 , respectively 
(in the MOND regime) is 

<j>(r) = (MGa ) 1 / 2 Zn j r ~ ri j . (31) 

(vii) More generally, all the physics of a system of gravitating point 'charges', q i} is encapsuled 
in the iV-mass energy functions E(r 1 , ...,r N ), since the forces are derived from them through F t = 
—dE/dri (The force on a given mass measures the change in the energy under a rigid translation 
of the mass.) 9 The above results give us the analytic expressions for E for the general two-body 
system: 

E(r 1 ,r 2 ) = ln\r 1 -r 2 \ f> ™, (32) 
and for a three-body system with vanishing total 'charge': 

E(r u r 2 ,r 3 )=ln[rf^ 3 rH% (33) 

where r id = \r t — r^| [Relation (32) is, in fact, a special case of this with a renormalized energy 
for the limit where mass 3 is sent to infinity (r 3 — > oo)]. We do not know the general form of E. 
For systems with vanishing total 'charge' we can show that E(r 1 ,...,r N ) is determined up to an 
unknown function of the N(N — 3)/2 conformally invariant variables u ijkl = r^r^jr^r^: 

E(r u r N ) = E (u ijkl ) + ^ b%ln r ij: 

l<i<j<N 

where (for TV > 2) 



(35) 

with v(q) = (2/3)(Ga ) 1 /2| g |3/2_ For iV = 3 we have no u ijkl variables, so Eq is a constant and 
Pij = bij. Both E and Eq have a suppressed dependence on the charges. 10 The second term in 



9 The gravitational energy of an isolated system with non-vanishing total 'charge' is infinite (it is finite for a 
vanishing total 'charge'), but its difference for two systems with the same 'charge' are finite, and only these interest 
us. 

10 In D dimensions u(q) <x q ^ - 1 '' . In two dimensions we have the exact result S(r 1; r N ) <x ^2 iK qiqjln r i} . 



(34) 



b N = — 

13 {N 



(iv-i; 



JV 
m= 1 
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eq.(34) carries the anomalous transformation properties of the energy under dilatations and special 
conformal transformations, while Eq is truly invariant. 11 For example, under dilatations 

N 

£(Ar!, Xr N ) = E(r u r N ) -InA^ v{q m ). (36) 

m— 1 

Again, by sending one of the charges to infinity (call it q N = —Q), and renormalizing the energy by 
subtracting the constant ln(r N ) b? N = —2u(Q)ln(r N ), which goes to infinity in the limit, we 

can define the energy for a system of charges with a finite total charge Q; for example, a system with 
only matter masses (good only for comparing systems with this total charge). The sum in eq.(34) 
then goes only up to N — 1, and the u variables with index ./V are written, e.g., u ijkN = r i} jr ik \ so 
they are still invariant under dilatations (but not inversions); so Eq is invariant under dilatations 
of ri, r N _ l . The energy for the remaining N — 1 masses transforms under dilatations as: 

E(Xr 1 ,...,Xr N _ 1 ) = Efa, r N _,) + Aln A, (37) 

where 

AT-l 

l<i<j<N-l m=l 

this is equivalent to eq.(27) (seen by taking the A derivative at A = 1). 

(viii) In general, the forces in two systems that are related by a conformal transformation of the 
charges are simply related (see Milgrom 1997). For example, consider a point charge q at a distance 
r < a from the center of a spherical shell of radius a, uniformly charged by —q. On dimensional 
grounds, the force on the point mass can be written as f(r) = —2v(q)ra~ 2 s{r/a). Then, from 
conformal invariance, the force on q when it is at R > a is f(R) = 2v(q)R~ 1 [\ + (a/ R) 2 s(a/ R)\ (the 
function s is not known). 

Under conformal transformations, equipotential surfaces go to equipotential surfaces of the 
new configuration, and field lines go to field lines. Also, spheres (circles) go to spheres (circles), 
including planes (straight lines), which are spheres of infinite radius. For an arbitrary system of 
charges, q,, with ^ ^ = 0, lying on a sphere of radius a, one can show 12 that the radial component 
of the force on q % is u{q i )/a (pointing outwards), and everywhere on the sphere, outside the charges, 
the field lines are tangent to the sphere (no radial force on test charges). 



The energy of a zero-total-charge system is invariant to all conformal transformations of the charge distribution. 
However, here we speak of transforming the point charges rigidly to their new positions, without affecting the 
transformation on their internal structure. The anomalous transformation properties of E(ri, r N ) result from this 
and the fact that the energy of a finite-charge system does change under dilatations, r — > Ar, by u(q)ln\. 

12 This is done by transforming the sphere into a plane where the forces must all lie from symmetry, and using the 
transformation law for forces from Milgrom (1997). 
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3.3. Well mixed configurations 

In subsection 3.2, I considered mainly systems of point masses for which matter and TM are 
well separated. It is worth noting that interesting situations may occur when matter and TM are 
well mixed in the sense that p « p. 

Look, for example, at the case p = p: In this case <fi = 0, and <fi is twice the Newtonian potential 
for each mass separately, call it 4> N . So the two sectors see the same potential, 4> = cf> = 2(<fi N . 
Given also the same initial velocities, such mass configurations will retain their density equality with 
time, with each mass type developing according to Newtonian dynamics with an effective Newton 
constant G//3. This echoes the situation in symmetric cosmology, where we have p M „ = g^, with 
each satisfying Freedmann's equations [possibly with a cosmological constant ~ A4(0)og] and with 
an effective gravitational constant G//3. Such local configurations are, however, unstable, and a 
small departure from equal densities will lead to eventual segregation of the two mass types. For 
example, in the presence of a small matter, or TM, mass p* (or with a small departure from p = p) a 
nonzero <f> potential is created which is added to <p and <p with opposite signs, causing the separation 
of p from p in a way that increases the separating force even further. Thus initially, at least, even 
if the surface densities of p and p separately are large compared with a /G, so that each mass type 
creates a high-acceleration Newtonian field, the resulting evolution of the difference 5 = p — p will 
be governed by MOND, and is nonlinear. 



3.4. Relativistic point mass solutions 

Allowing for TM in BIMOND permits a two-parameter (M, M) family of (relativistic) spher- 
ically symmetric, static, vacuum solutions on a double Minkowski background (neglecting 'cosmo- 
logical constant' effects). These would describe, e.g., matter-TM black holes. For M = M we 
have the solutions <? M „ = g^, with both metrics being of the Schwarzschild form for M, with G//3 
as gravitational constant. At the other end, for M = we can have pure-matter black holes, for 
instance. Because of the cosmological coincidence a cHq/2tt, for all sub-Universe systems the 
horizon acceleration is much larger than o . This means that MOND effects enter only far beyond 
the horizon, deep in the NR regime. For /3 = 1 we know that in the high-acceleration regime 
BIMOND describes two separate Einstein theories for matter and TM [with a cosmological con- 
stant ~ a„A4(oo), which I neglect here]. In this case, each of the metrics is nearly a Schwarzschild 
one corresponding to its own mass, within many Schwarzschild radii. At large radii the metrics 
approach their NR expressions as per eqs.(4)(23-24): at first <p « —MG/r, ip rs —MG/r, and only 
at much larger radii (j) w -0 w 4> k, [(M - M)Ga ] 1/2 ln r (for M > M). 
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3.5. Asymptotic field of equal-mass systems 

For an isolated system with equal total masses of matter and TM, M = M, the asymptotic 
behavior of the MOND potential is not described by the usual logarithmic dependence on the 
radius. This can be easily seen by applying the Gauss theorem to eq.(17), which implies a vanishing 
coefficient to a logarithmic term. Instead, the generic asymptotic potential is of the form 

0(r) -> (MGa„) 1/2 R- ^, (39) 

where R is some radius vector characteristic of the mass distribution in the system. This is the 
vacuum solution of eq.(17) with the slowest decrease with radius. It is gotten by a conformal 
transformation (inversion at the origin) from the configuration of a constant acceleration field 
4> oc R • r, which is clearly a vacuum solution, and is thus itself a vacuum solution by virtue of the 
conformal invariance of eq.(17) (Milgrom 1997). 13 Not surprisingly, the potential in eq.(39) has 
the form of a dipolar electrostatic field in two dimensions: D • r/r 2 . However, here the asymptotic 
"dipole strength" is not related to the dipole of the mass distribution; in fact, I do not know 
how to express R in terms of the mass distribution. For a matter-TM point-mass pair, the exact 
expression, eq.(31), implies R = r 2 — r 1 . More generally, we know that R does not scale with the 
mass, and scales with the system size: when p(r) — > A _3 p(r/A), we have R — > AR. 

The magnitude of the asymptotic field in eq.(39), \V<p\ = (MGa ) 1 / 2 |R|r -2 , depends only on 
r (not on the direction) and decreases with radius like the Newtonian field of the system \Vct>\ = 
2MGr~ 2 . Their ratio is \V(j>\/\V(f)\ = 2i? M /|R|, where R M is the MOND radius of the system. 
So, if the mass distribution is such that |R| 3> R M , the MOND contribution dominates the field 
asymptotically as well. 

There are configurations for which |R| is much smaller than the characteristic size of the 
system. In particular, there are systems for which R = 0, for example due to symmetry. For these, 
the MOND potential decreases faster asymptotically, and the Newtonian field dominates there. I 
have not identified the next slowest decreasing vacuum solution of eq.(17). Because of the conformal 
invariance of the problem, finding such solutions is equivalent to finding the behavior of vacuum 
solutions near the symmetry point of a symmetric mass configuration. Take, for instance, a system 
with two matter masses, M, at opposite corners of a square, centered at the origin, and two TM 
masses M = M at the other two corners. Inversion about a sphere centered at the origin, and 
containing the corners, leaves the mass configuration invariant, but interchanges the origin with 
infinity. This implies that the behaviors of the potential near the origin and at infinity are gotten 
from each other by a simple inversion. Spherically symmetric configurations with M = M have a 
vanishing MOND potential (f> outside the mass, with only <p contributing. 



The logarithmic vacuum solution transforms to itself. 



-16- 



4. Observational consequences 

No BIMOND effects are felt in the cosmos, in the picture in which on cosmological scales 
% v = as long as matter is homogeneous (except perhaps those of a cosmological constant). 
BIMOND effects appear when inhomogeneities develop-the beginning of structure formation-when 
the presumably uncorrelated perturbations induced in the two sectors would have caused local 
departures from equality of the metrics. The initial stages of the growth of perturbations on the 
background of the expanding Universe would have been affected already by some of the peculiar 
matter-TM gravitational interaction discussed here, but these require using the relativistic version 
of the theory. Subsequent, nonlinear development of inhomogeneities occur already in the NR 
regime, for which we can directly use the NR theories, and the results of this paper. Matter-TM 
interactions can also lead to observable effects in the present day Universe. Here I discuss very 
succinctly some possible consequences of the existence of TM and its BIMOND interaction with 
matter. 

4-0.1. Structure formation 

Even if the underlying theory itself is symmetric in the two sectors, it is not clear that the 
attributes of their cosmic matter contents are exactly the same. For example, we do not yet know 
what engendered the baryon asymmetry in matter. It is possible, then, that the small initial 
asymmetry was different in the two sectors, which could result in different baryon densities today. 
Also, even if there occurred similar inflation episodes in the two sectors, they could have emerged 
from them with different perturbation amplitudes. The effects of such possible asymmetries, and 
others, on structure formation should be studied. 

The presence of TM and its peculiar interactions with matter, discussed here, would have 
had major impact on the process of structure formation. As already mentioned in the introduc- 
tion, in a symmetric universe, with g^ v = g^ v on large scales, small fluctuations are governed by 
nonlinear MOND dynamics, even at early times when the universe at large is characterized by 
high accelerations, cH S> a . What dictates that we are in the MOND regime is the smallness of 
the perturbations, leading to small departures from metric equality between the two sectors: The 
argument of the interaction function Ai, or its derivative M! appearing in the field equations, is 
the difference in the connections of the two metrics in the relativistic regime, or the difference in 
the gradients of the matter and TM potentials in the NR regime, in units of a . The growth of 
perturbations in this picture is thus nonlinear from the start, since the MOND potential is nonlin- 
ear in the overdensity. To apply the present results to the problem would thus require numerical 
simulations. 

Simulations of structure formation with NR MOND dynamics has been considered in several 
studies; e.g., by Sanders (2001), Nusser (2002), Stachniewicz & Kutschera (2002), Knebe (2005), 
and Llinares, Knebe, & Zhao (2008). It should be easy to extend such simulation to reckon with 
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TM, by including an initial TM distribution with its own seed fluctuations, and employing the field 
equations (14)(16), or alternatively the more wieldy eqs.(22)(16). 

Our results here are not directly applicable to the early stages of the growth of perturbations, 
because they assume that the system at hand is much smaller than the characteristic cosmological 
curvature radius, and they ignore the underlying cosmic expansion, relying on perturbations around 
Minkowski metrics, not those describing cosmology. However, these results would still apply at later 
times, and might be indicative of what happens qualitatively even at early times. 

One important effect that can be anticipated is the segregation of matter and TM due to their 
mutual repulsion. This would presumably lead to formation of interweaving cosmic webs for matter 
and TM; i.e., complementary networks of nodes, filaments, and voids, with the high-density nodes 
of one sector residing in the voids of the other, and vice- versa, with filaments avoiding each other. 

Some of the main questions that can be answered by numerical simulations are: Do matter and 
TM indeed segregate efficiently and what does the present-day configuration look like? Do they 
indeed form interweaving, mutually avoiding cosmic webs? Is segregation practically complete, or 
do we still find galaxy size or larger TM bodies lingering in the neighborhood of matter structures? 
Does this occur often enough for direct gravitational effects of TM on matter to be a common 
phenomenon? To what extent voids are more pronounced (empty) compared with simulations with 
matter only? What are the effects on peculiar motions of large-scale structures? 

Regarding the question of void structure, in particular, it has been suggested that various 
observed aspects of voids pose problems for the LCDM paradigm (e.g., Peebles 2001, Tully 2007, 
Tully & al. 2008, Tikhonov k Klipin 2009, Tikhonov & al. 2009, and see van de Weygaert & 
Platen 2009, and Peebles & Nusser 2010 for recent reviews). Even without TM, MOND is expected 
to produce more pronounced voids than DM for the following reason: A void acts like a region of 
negative mass. If the ambient density is p a , then a void can have at most the effective density — p a . 
In MOND this affective density can be much higher if the characteristic accelerations involved are 
small enough. TM can help even more because of the added repulsion by TM concentrations, which 
would arguably reside in matter voids. 

Simulations should also study the effects of various possible departures from symmetry between 
the two sectors, such as disparate baryon densities, different amplitudes of the initial fluctuations, 
etc. One should also investigate the dependence of the outcome on the parameter £, and also 
consider structure formation in asymmetric theories with a ^ f3. 

4-0.2. Gravitational lensing 

The best prospects for detecting or constraining TM, if it has formed structures like those of 
matter, and if it has segregated efficiently from matter, seem to be via weak gravitational lensing. 

Expressions (4) for the two NR metrics imply that we can calculate lensing in the standard 
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way, as is done in GR, using the MOND potentials instead of the Newtonian potentials. 14 

We are interested here in the lensing effects of an isolated TM body, p, on matter photons. 
These are dictated by the matter MOND potential 4> created by p. Use, as an example, the field 
equations in the form of eq.(8), with a = (3 = 1, and p = 0. We then have 

= -4irG Pp , = 47rG(p + p p ), (40) 

where p p = — V- (A^'V0*). In other words, we can calculate (j) as the Newtonian potential produced 
by minus the phantom density produced by p alone. This phantom density is the same as we would 
calculate for a matter distribution p = p. Thus lensing of matter photons by a TM body are 
different in two major ways from lensing by a matter body of the same mass distribution: First, 
photons (or any matter test particle) do not 'see' the TM itself, only the fictitious, MOND, phantom 
matter that it produces. Second, photons sense this phantom matter with an opposite sign: a TM 
lens thus acts as a diverging lens. 

The first fact above implies that photons are oblivious to the existence of the TM in regions 
of high acceleration: For a spherical body with radially decreasing acceleration, no force is felt in 
the Newtonian regime. For example, a point TM mass M, is not felt by matter roughly within its 
MOND radius R M = {MG/a Q ) 1 / 2 . 

Another fact to recall is that the thin-lens approximation does not apply in MOND (Mortlock 
and Turner 2001, Milgrom 2002b, Milgrom &: Sanders 2008) and structure along the line of sight 
has to be reckoned with: What enters lensing is still the integrated surface (column) density of the 
phantom matter. This, however, does not depend only on the surface density distribution of the 
baryonic matter that produces it; it also depends crucially on the baryon distribution along the 
lines of sight. For simplicity's sake, I consider here only a spherical lens. In this case, the ratio of 
the Einstein radius of a lens to its MOND radius is 



with d l3 , d t , d s the lens-source, lens, and source angular-diameter distances, respectively. Consid- 
ering the well known proximity a pa cH /2ir, this ratio is 



R M \ird s D H 



(42) 



where D H = c/Hq is the Hubble distance. Thus, for lenses nearer than the Hubble distance, the 
Einstein radius is within the MOND radius. But no bending occurs within the MOND radius, since 
M itself is not felt, and there is no PM there. Also, for a spherical lens, the PM is by definition 



14 This is not the case for BIMOND theories that hinge on scalar arguments of the bimetric interaction other than 
that given in eq.(3). 
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always at or outside its MOND radius and so definitely outside its Einstein radius (the surface 
density of PM is always subcritical for a spherical body-see below). We thus do not expect strong 
lensing effects of TM on matter photons. 15 

Since only the phantom matter affects lensing of matter photons by TM, we should recall some 
of its important properties. I mentioned already that it is to be found only roughly beyond the 
MOND radius. Another pertinent attribute is that its local surface density (for a spherical system) 
is bounded by some universal value: S < So ~ a /2irG. This is related to the observation of 
Brada & Milgrom (1999) that the acceleration produced by the MOND PM can never much exceed 
a . More directly, I showed in Milgrom (2009a) that the MOND prediction for the central surface 
density of phantom haloes in spheroidal systems is < So- This fact is also supported by observations 
of galaxies (Donato & al. 2009, Gentile & al. 2009). The exact value of the maximum as deduced 
from MOND depends somewhat on the MOND formulation and on the choice of extrapolating 
function. In comparison with the lensing critical surface density of the lens E c = c 2 d s /4TrGd ls d l we 
have, using again the numerical proximity of a ps cHq/2tt, 

— = — l ±J__ (43) 

So PM lensing, and as a result TM lensing of matter photons, is always sub-critical for spherical 
lenses, and very much so for low redshift lenses. 

A possible, interesting signature of lensing by the PM of a TM body is the appearance of a 
ring, or a projected shell of PM (negative mass in our case) if the TM body is contained within its 
MOND radius, as discussed in Milgrom & Sanders (2008). With lensing by a matter body (such as 
a galaxy cluster) the weak ring feature is generally masked by the dominant surface density of the 
baryonic body. But in the case of lensing by TM, the TM baryons do not contribute to lensing, 
and thus do not mask the feature if it is produced. 

We should try then to identify the signature of possible TM bodies, such as TM galaxy clusters, 
in weak lensing surveys. Several factors may turn out to make such identification difficult: In the 
first place, we can, of course, not rely on guidance, or substantiating evidence, from matter effects- 
as can be done in standard weak lensing-such as x-ray emission, Sunyaev-Zel'dovich effect, or 
a direct view of a galaxy concentration. Second, as we saw, no high surface density effects are 
expected, since only the rather law-surface-density fictitious PM acts as a lens. Third, the signal 
should evince the effects of a negative mass concentration, which may require specialized algorithms 
to discover. Fourth, matter structures along the line of sight will act to cancel the lensing effects 
of TM bodies, or at least distort and confuse them. 



15 With a TM lens that is elongated along the line of sight we can have strong lensing. For example, consider N 
equal point TM masses, M, along the line of sight, separated from each other by more then their individual MOND 
radius R M - The projected PM is then N times the individual contribution, so the Einstein radius of the system is 
N 1 / 2 times that of a single mass. For large enough N, this is larger than R M , and we can then get strong lensing 
effects. Also, with /? / 1 there will be attractive or repulsive strong lensing effects. 
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Results of large-scale-structure simulations with TM will pinpoint the expectations for weak 
lensing by TM, and help design search strategies. 

4.0.3. Other effects 

If matter-TM segregation is not complete, and galaxy-size TM bodies are still lurking in the 
neighborhood of matter galaxies, they may produce distortions that are visible as warps, lopsid- 
edness, or ellipticity. In particular, in theories with /3 < 1, matter and TM attract each other 
in the high acceleration regime; so, there might be TM trapped in high acceleration regimes of 
matter territory. One may speculate in such a case that small amounts of TM trapped in the (high 
acceleration) cores of galaxy clusters might be responsible for the observed mass discrepancy there 
(see 5.3 below). 

If structure formation simulations confirm that indeed matter and TM form interleaving mu- 
tually avoiding webs, and if there is matter-TM symmetry in the cosmos, we can estimate the 
occurrence of TM galaxies in matter territory, by seeing how often we find matter galaxies in voids, 
which are presumably TM territories. 

Clearly, there should also be important effects of the existence of TM on the appearance of 
the CMB fluctuations. 

5. Asymmetric theories 

In theories with a 7^ /3 we have different gravitational dynamics in the two sectors. In partic- 
ular, in such theories the amount of TM in the universe, and its general properties (distribution, 
etc.) can be rather different from those of matter. Such theories are, clearly, worth investigating, 
even if they involve treating less amenable configurations. Here I discuss briefly the NR limit of 
several examples of such theories. 

5.1. A theory with j3 = 1 and a 2> 1 

In this example we take the extreme limit a — > 00 while (3 = 1. We can then write the field 
equations for the MOND potentials, eqs. (14-16), as 

A(C<^) = 4nGp, V • {M'[(V0/a o ) 2 ]V^} = AkG( P - p), (44) 

4> = (4> + 4>, 4> = (4>- (45) 

Since now £ — > 0, we have TW'(oo) = 1. We find then that (,(f) = (fr N is the Newtonian potential of 
p, while 4> is a solution of the nonlinear Poisson equation V • [p(\V(j)\/a )V4>] = Ai:G{p — p), with p 
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the standard MOND interpolating function, and <j> = <f) + 4> N , (f> = <j> N . Dynamics in the twin sector 
is thus fully Newtonian and oblivious to the matter sector. In contradistinction, the gravitational 
potential in the matter sector is governed by the sum of the MOND potential produced via the 
nonlinear Poisson equation by p — p, with the Newtonian potential of twin matter. 

5.2. A theory with a + /3 = 

As another example, take the case a + j3 = 0, which leads to the quasidinear formulation of 
MOND (QUMOND) discussed at length in Milgrom (2009c). The field equations are then 

= 47TG/TV + /r x V • (M'V(f>*) = 47tG( P + p p ), 

A<f> = 47rG/rV + /r x V • (M'V</>*) = 4irG(p + p p ), (46) 

where, (f>* = <j> — <j>, and 

Pv = (4^G)- 1 /3- 1 V • (M'V<f>*) + (r 1 ~ 

p p = (4^G)- 1 /3" 1 V • (M'V<P*) + {p' 1 ~ l)/5- (47) 
Taking the difference of these equations gives 

A(f>* = 47rG/3~ 1 (p-p). (48) 

One first solve the Poison eq.(48) for <j>*, and then another Poisson equation for <\> or <p. Accelerations 
of test particles are given by a = — V</>, a = — V<f>. 

We see that, again, the dynamics within each sector separately are not the same. If we compare 
the NR potentials for two configurations p = p R , p = 0, and p = p R , p = 0, we see that cf>* has an 
opposite sign for the two configurations, and hence p p / p p . 

The conserved momentum is now P = f (pv — pv). If we still define the force on a matter 
subsystem in the volume v as F = — f v pV(f>d 3 r, and that on TM as F = — pV(j)d 3 r, then for a 
closed system made of matter and TM it is F — F that vanishes, not their sum. If they do not vanish 
separately, such an isolated system will self accelerate (still preserving momentum). This seems 
paradoxical if applied to an isolated system, but in the context of a universe filled with matter 
and TM in equal amount it does not necessarily lead to unacceptable behavior since no large scale 
directed accelerations can occur. It remains to be checked by numerical simulations whether this 
can lead to inconsistencies (conceptual or observational) in the expected behavior of matter. 

Specialize further to j3 = 1, which is particularly transparent: In this case eq.(48) implies 
that (f>* is the Newtonian potential of the system with TM contributing as having negative (active) 
gravitational mass. In eq.(47) we have p p = p p = (47rG') _1 V- (A4'V(j)*). To get a Newtonian limit in 
the matter sector we have to have M'{z) — s> for z — > oo [z = (V^*/a () ) 2 ]. This also gives standard 
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Newtonian dynamics in the TM sector. It also means that matter and TM do not interact in the 
deep Newtonian regime. In this case the total force on matter and on TM for a closed system 
vanish separately. 

To get MOND dynamics in the matter sector we have to have M'{z) ~ z~ 1/4 for z < 1. 
Thus effectively, it can be said that we can use a Newtonian calculation with each type of mass 
'seeing' the density produced by its own type plus the phantom mass that is common to both 
sectors. Because p enters the source of 0* with an opposite sign one may say that the phantom 
mass produced by TM alone is repulsive to both types of matter, while that produced by matter 
is attractive to both. Thus, in the TM sector, a point mass M produces a Newtonian, attractive 
force on a TM test particle within its MOND radius, R M , but instead of the MOND enhancement 
of gravitational attraction in the matter sector, in the TM sector MOND effects cause this force to 
weaken the attraction, and turn it into repulsion at large distances-a manifestation of the disparate 
behavior in the two sectors. 



5.3. A theory with < (3 < 1 and a > 1 

Here I consider an example with < j3 < 1, and to accentuate matters, with f3 <C 1. In such 
theories matter and TM attract in the high acceleration limit but still repel in the deep-MOND 
limit. Since we have to have ( < 1 we I take a > 1, or, for concreteness' sake take a = 1. 

From eq.(16) we can now write the potentials for the two sectors as 

<p = 4> N + /r 1 ^ + j>, $ = cf> N + r l 4> N - ^, (49) 

where I used the fact that A = fi/a <C 1, and <p is a solution of the second of eq.(14), i.e., 
V • {Ai 1 '[(V (ft I a ) 2 ]V '<fi} = 4irG(p — p). The gravitational potential that matter sees is thus the 
sum of three terms: its own Newtonian potential, the Newtonian potential of the TM enhanced 
by the large factor /3 _1 , and the "MOND" potential 4>. Note that because here ( « 1, we have 
M'(oo) » A" 1 > 1, so for a -> 0, (f> -»• A(0 JV - <f> N ). 

We see that even a small admixture of TM in an otherwise pure matter object can make an 
important contribution to the matter potential, because its Newtonian potential is enhanced by the 
factor /3" 1 . For example, in such a theory, a small amount of TM trapped in the cores of (matter) 
galaxy clusters could explain away the mass discrepancies observed there. Such amounts of TM 
need not affect much the behavior in the MOND regime of the clusters, since TM enters the source 
of the second of eq.(14) with the same weight as matter. 

It remains to be seen whether such asymmetric theories can be made consistent with cosmology. 
Recall that it is G//3 that plays the role of the Newton constant in the BIMOND cosmologies I 
have considered so far. This can be interpreted as an apparent enhancement of the matter density 
by a factor For example, in the Friedmann equations a baryon density p b appears as p b //3, 

which would be interpreted as baryons plus DM of density /? b (/3 — 1 — 1). 
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6. Summary and discussion 

I have studied some aspects of the dynamics of matter an the putative TM that may be 
present in the context of BIMOND. This is done for NR bodies of either type. In fully symmetric 
theories, which seem preferable on various grounds, we get MOND dynamics within each sector. 
The interaction between matter and TM is, however, nonstandard even compared with MOND: In 
the deep-MOND regime matter and TM bodies repel each other with MOND-like forces (decreasing 
as inverse distance). In the Newtonian, high-acceleration regime the force depends on the parameter 
/3. For the fiducial value (3 = 1, the matter-TM interaction vanishes. For /3 < 1 there is Newtonian- 
like attraction (decreasing as inverse squared distance) with strength /5 _1 — 1, and for /3 > 1 there 
is similar repulsion. 

I have also considered briefly possible effects of the presence of TM with such properties on 
structure formation and its lensing properties. To assess more reliably such effects we have to call 
upon numerical simulations. 

It needs to be stressed, finally, that the TM does not, indeed cannot, play the full role of dark 
matter in galactic systems (although, as we saw, it may produce some effects attributed to dark 
matter): There isn't enough of it; it probably shies matter galactic systems; and, in all probability, 
would, anyhow, decrease gravitational attraction between matter bodies if it comes between them. 
It is still the MOND departure from standard gravity that is responsible for the observed mass 
discrepancy. 
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